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Abstract 

We analyze the dynamics of warm inflation models with general viscous effects. We characterise 
"^j." ■ the situations yielding asymptotic scaling behaviour. 
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I. INTRODUCTION 



Scenarios leading the universe to a moderate temperature state at the end of the inflation 
avoiding the need of a separate stage of reheating [l[ |^ have been proposed by several 
authors. It was advocated that this can be accomplished by coupling the inflaton to the 
radiation field in such a way that the decrease in the energy density of the latter during 
inflation is somewhat compensated by the decay of the inflaton into radiation |3|. This 
kind of mechanism, first proposed by Berera was termed "warm inflation" since the 
radiation temperature never drops dramatically. In addition, it has two further advantages: 
(i) the slow-roll condition <fi 2 <C V{4>) can be fulfilled for steeper potentials, and (ii) the 
density perturbations generated by thermal fluctuations may be larger than those of quantum 
origin [5|. 

One natural further step is to assume that the inflaton field decays not just into radiation 
but into massive particles as well. The mixture of massless and massive particles can be 
described by an overall fluid with equation of state p = (7 — l)p, where the adiabatic index 
7 is bounded by 1 < 7 < 2. It is also natural to expect that this fluid has a negative 
dissipative pressure, II, that quantifies the departure from thermodynamical equilibrium. 
This viscous pressure arises on rather general grounds from interparticle interactions 

Our purpose in this work is to analyse the dynamics of models that generalize the usual 
warm inflationary scenario by introducing these novel elements, namely the decay of the 
scalar field into a fluid of adiabatic index 7 rather than just radiation, and specially the dis- 
sipative pressure of this fluid. We will not dwelve into the difficult question of the quantum, 
nonequilibrium thermodynamical problem underlying warm inflation [7L but rather take a 
phenomenological approach similar to that considered in several works jD, 0] . We allow for 
various forms of the rate of decay of the scalar field, as well as for various expressions of the 
dissipative pressure, and resort to the qualitative analysis of the corresponding autonomous 
system of differential equations using the approach developed in This enables us to 

consider arbitrary scalar field potentials. 

A central question we address is the existence of scaling solutions 0,[l(J U|, i-e., solutions 



where the ratio of energies involving the matter fluid and scalar field keep a constant ratio. 
In the standard treatment of warm inflation the existence of a "quasi-static" regime for the 
energy transfer between the scalar field and the radiation is assumed for the computation of 
the temperature of matter at the final stages of warm inflation (Hf . It is however unclear why 
there should be such a quasi-static stage in the conventional model of warm inflation, and 
thus whether such an assumption is fully legitimate. If there were to be a scaling regime at 
the end of warm inflation such a quasi-static regime would necessarily exist. Moreover, this 
class of solutions finds itself a justification in the kinetic analysis of interacting fluids [l2^ . 
A particular illustration of this can be found in 

In this work we show that the existence of a late time scaling behaviour depends not only 
on the asymptotic form of the inflaton potential [? ], but also on having a time-dependent 
rate of decay for the scalar field. 



II. THE DYNAMICAL SYSTEM 

We consider a spatially flat Friedmann-Robertson- Walker Universe filled with a self- 
interacting scalar field and a perfect fluid consisting in a mixture of matter and radiation, 
such that the former decays into the latter at some rate T and that the latter exhibits a 
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dissipative pressure II. On the other hand, we neglect radiative corrections to the inflaton 
potential. The corresponding system of equations reads 

3H 2 = p+^- + V{4>) (1) 

H = -\ (> + 7P + n) (2) 

4> = -(3H + r)4>-v'(4>) (3) 
p = -3^( 7 p + n) + r0 2 , (4) 

where we have used units such that 8irG — c 2 — 1. The first two are Einstein's equations, 
the third describes the decay of the inflaton, and the fourth one is the energy balance for 
the matter fluid. As usual H = a/a denotes the Hubble factor. 

To cast the corresponding autonomous system of differential equations it is expedient to 
introduce the set of normalized, dimensionless variables 

x 2 = ^- V * = YM_ r = JL (5) 

QH 2 ' y 3H 2 ' 3H' { ' 

as well as \ = U/3H 2 , along with the new time variable N = In a. We assume, quite 
generally, = T((p)H s and IT = — 3^p a H 13 , where 5, £, a and /3 are constants and 
moreover 2a + (3 — 2 = on dimensional grounds. The latter condition on the parameters 
a and j3 implies that the dimensionless variable x becomes x = ~3 a ( (p/3H 2 ) a , and hence 
reduces to x — ~3 a ( (1 — x 2 — y 2 ) a , since p/3H 2 = 1 — x 2 — y 2 . Therefore the dynamical 
system is 4-dimensional and reads 

x > = x (Q-3(i + r ))- W {<t>)y 2 , (6) 
V = (Q + W(<P)x)y, (7) 



v^l^l x + Q(l-6) 



(8) 



r 

0' = V6x , (9) 
where a prime indicates differetentiation with respect to N, and where we have used the 
definitions W(<f>) = ^/| and 

Q = l [2x 2 + 1 (l-x 2 -y 2 )~3 a ai-x 2 -y 2 r] . (10) 



Equation @ was first considered in jipj] and is crucial for the consideration of general 
potentials V(<p) besides the particular case of the exponential potential. The function Q 
defined by Eq. (fTUj) is related to the deceleration parameter q = —da/ a 2 by q = Q — 1. 

III. QUALITATIVE ANALYSIS 

A. The special case oc H and II = 

Here we take = 3rH, where r is a positive constant. The dynamical system reduces 
to the three equations 

x' = x [Q_3(l + r)] -W(<P)y 2 , (11) 
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V = [Q + W(</>)x]y, (12) 
(f)' = V6x , (13) 

where Q is now simply given by Q = | [2a; 2 + 7 (1 — x 2 — y 2 )]. 

We see that these equations are analogous to those of the r = case with a different 
coefficient on the linear term in x of Eq. (|llj). Thus the basic qualitative dynamical features 
remain the same as those found for that model 0, 0, The decay of the scalar field 
introduces though two major consequences. 

The origin x = 0, y = is again a fixed point associated with the vanishing of the 
scalar field's energy and, hence, corresponding to matter domination, but the interaction 
given by a non-vanishing r has the relevant effect (already found in the constant T model) 
that the stability of the minima is reinforced and that the maxima become less unstable. 
The other major effect of the interaction arises when we look for fixed points at <p — > 00. 
Now, we find that there are always attracting scaling solutions for potentials that have an 
asymptotic exponential behaviour, that is, for potentials for which W const when <p — > 
00 (lfj. Moreover, this happens independently of the steepness of the late time exponential 
behaviour, a remarkable effect of the present model. 

Assuming the potential to be asymptotically given by V oc e~ x< ^, these solutions are 
a(t) oc t A , (p — 0o = lnt ±2 / A , where A is given in implicit form by 



37U--)U-T^(l + r) -^ = 0. (14) 



A linear expansion in r in the neighborhood of the scaling solutions of the r = case yields 



A = — 

3 7 



1 + 



2_ 

A 2 



2 2_ 

3 7 A 2 



(15) 



So the decays have the effect of increasing (resp. decreasing) the scale factor rate of expansion 
with regard to the r = case if A 2 > 37 (resp. A 2 < 37). In particular we can see that the 
scaling behaviour can be inflationary, for cases where this would not happen in the absence 
of decays. For instance, taking 7 = 4/3 and A 2 > 4 the condition for the scaling solution to 
be inflationary is 1 + r > A 2 /4 > 1. Thus, in this model, the solutions yield endless power- 
law inflation even for a modest scalar field decay, provided that the asymptotic behaviour 
of the potential is steep enough, i.e., A 2 > 37 (> 4 in the present example). Besides the 
scaling solutions, there are also fixed points corresponding to de Sitter behaviour x = 0, 
y = 1, whenever the scalar field potential exhibits an asymptotic, non- vanishing constant 
value. However, when the potential is asymptotically exponential, there are no fixed points 
on the boundary x 2 + y 2 = 1 at — > 00, labelled 0oo, in contrast to what happens in the 
r = case. 

From a thermodynamical viewpoint the former scaling solutions are particularly interest- 
ing. In a universe with two components it can be shown [l2| that the temperature of each 
of the components satisfies the equation 

Tj _ a ( Tj \ dpi rijSj 

T % a V 3HJ d Pl dpi/dTi ' 1 ' 

where (i = 1,2), rij denotes the number density of particles of the z-species, Tj their rate 
of change, and Tj the temperature of this component. In the important case of particle 
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production happens with constant entropy per particle, = 0, we also have (pi +pi)Ti = 
— (p2 + P2)^2- Thus, taking the first component to be the matter/radiation fluid, and the 
second to be the inflaton scalar field, we have 

(p + p) . . 

J- 4> = ^7 J- m/r OC 1 m/r- {± ' ) 

As (p+p)/0 2 = 7(1 — x 2 — y 2 )/x 2 is a constant at the scaling solutions, then = 3rH, with 
r a constant, implies T m / r = 3a H, where a is another constant that depends both on r and 
on the location of the scaling solution. This yields a temperature of the matter /radiation 
component evolving as a power-law T oc a~( ry ~ 1 ^ 1 ~ a > . Thus for a close to 1, the temperature 
of the matter/radiation remains quasi-static, whereas for o > 1 (resp. a < 1) it increases 
(resp. decreases) . Provided that we guarantee enough inflation, r need not be very large 
(contrary to what is usually assumed to facilitate the slow-rolling). 



B. The general case 

In what follows we return to the general system (JHJl-Q- It seems reasonable to further 
assume that < a < 1 so that (3 > which amounts to having a bulk viscosity pressure 
whose importance diminishes with the expansion and with the dilution of matter. 

At finite values of we find a line of fixed points x = y = r = 0(<f)is arbitrary), and a 
line of fixed points characterized by x = 0, = 0o, (1 — U 2 ) l a — 3"C/7 an d any value of r. 
In the latter case O is the value of at an extremum of V((p), i.e., where W(4> ) = 0, and 
we require that 

3°C < 7 • (18) 

The linear stability analysis shows that the singular points x = y = r = corresponding 
to matter domination are unstable. Regarding the line of singular points with r ^ 0, linear 
stability analysis reveals that, besides the vanishing eigenvalue associated with r (A r = 0), 
the sability is once more determined by the nature of the extremum of V((f>). Indeed, one 
eigenvalue is always negative, and the eigenvalues along x and are given by 

A*,, = -Ml±l) ± i ^9(1 + r)* - AVQ y 2 W'(<M , (19) 

where is the equilibrium y-coordinate. We see from the latter equation that > 
requires that W(0*) < 0, that is a maximum at V(0*). Otherwise, in the case of a minimum 
of V(4>), we have either a stable node (when < W // (0 O ) < 9(1 + r) 2 /A\/Qy 2 ) or a stable 
sink (when W((f)o) > 9(1 +r) 2 /A^/6y 2 > 0). Decreasing y* induces an effective reduction of 
the steepness of the potential at the maxima and a greater stability of the minima. The net 
effect is that the system spends more time in the neighborhood of a fixed point associated 
with a maximum of the potential (alternatively, the minima become more stable). This is 
helpful for setting the conditions for slow-roll inflation, which translate as 3(1 + r) > y*, and 
do not require a large rate of decay of the scalar field. 

At 0oo, Eq. (jHJ) shows that we may have the usual fixed points corresponding to a non- 
vanishing, flat asymptotic asymptotic behaviour of the potential (late-time approach to a 
cosmological constant) if Q = 0, T'/T = and IT/II = simultaneously. However, we also 
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find asymptotic scaling behaviour in the case W((j)) approaches an exponential behaviour 
(W((j)oo) = —A, with constant A > 0), provided 

WW = (20) 

which amounts to having r(0oo) = v (V"(0 OO )) 2 , where v is a constant. We thus effectively 
have a oc H rate of decay, as in the previous subsection. Hence this emerges a necessary 
condition for scaling behaviour. In polar coordinates x = p cos 8, y = p sin#, these scaling 
solutions are given by 

cosfl* = , , y/u, (21) 
3(1 + r) 

where it = p 2 is a root of the quadratic equation 

(1 -u)(a- bu) - rbu - — (1 - u) a = , (22) 



and a = 2 and 6 = j^f^y- 

We see from Eqs. ()21j) and (j22j) that there is always one (and only one) scaling solution, 
provided the condition ()18j) holds. (Notice that this was precisely the condition that was 
required for the existence of fixed points at finite 0). We also find that the location of the 
singular points corresponding to the scaling behaviour is now closer to x 2 + y 2 = 0, having 
a smaller y^ value than in the models without bulk viscosity. Moreover, linear stability 
analysis reveals that under the conditions (fTHj) and a < 1 the scaling solutions are stable, 
i.e., they are attractors. These results mean that the late time contribution of the matter 
component is enhanced by the viscous pressure. This is a most convenient feature for the 
warm inflation scenario, since it further alleviates the depletion of matter during inflation 
and the subsequent need for reheating. 
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